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Abstract. We prove a single sum formula for the linearization co- 
efficients of the Bessel polynomials. In two special cases we show 
that our formula reduces indeed to Berg and Vignat's formulas in 
their proof of the positivity results about these coefficients (Con- 
structive Approximation, 27 (2008), 15-32). As a bonus we also 
obtain a generalization of an integral formula of Boros and Moll 
(J. Comput. Appl. Math. 106 (1999), 361-368). 
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1. Introduction 
The Bessel polynomials g„ of degree n are defined by 



fc=0 



where we use the Pochhammer symbol {z)n := z{z + 1) . . .{z + n — 1) 
for z G C and n G N. The first values are 

qo{u) = l, qi{'u) = l + u, q2{'u) = I + u + —. 

Using hypergeometric functions, we have qn{u) = iFi{—n;—2n;2u). 
They are normalized according to qn{0) = 1, and thus differ from the 
monic normalization 6n{u) in Grosswald's monograph [6]: 

(2n)' 

The polynomials 6n are sometimes called the reverse Bessel polynomials 
and yn{u) = u'^9n{^) the ordinary Bessel polynomials. These Bessel 
polynomials are, then, written as 



, . (2n)! „ ,1, ^ {n + k)\ , . , 

fe=0 ^ ^ 
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The so-called linearization coefficients /3^'"^^ of the Bessel polynomi- 
als [2] are defined by 



n+m 



qn{au)qm{{l - a)u) = ^ Pk'''''' {a)qk{u). 



For example, we have /3 



(3,5) 



j(n,m) / 

I- 
k=0 

for = 0, 1, 2 and 



(3) 



^(3,5) 


[a) 


= a\ 




^(3,5) 


[a) 


= 7a^(l 


- a) , 


^(3,5) 


[a) 




^ (84 - 126 + 126 - 84 + 36 - 


^(3,5) 


[a) 




- af (126 - 168 + 108 - 36 a + 5) , 


^(3,5) 


[a) 


143 2 
5 


(1 - a)^ (I2a2 - 9a + 2) , 


^(3,5) 


[a) 


= 143 





Recently Berg and Vignat [21 [3] have established some important results 
about the coefficients /3fc"'™''(a.) with applications. In particular, they 
proved [21 Theorem 2.1] that /3i"'°''(a) = a"" and, for < A; < n — 1, 



(n,0) 
re 

n\ (n-i-l)AA: 



/3r'°)(a)=a'=(l-a)-g- 

bfc/ j=o 



n + l\/?2 — /c — 1 
A; - j 



(l-ay. (4) 



Moreover, for n,m > 1, Berg and Vignat ^ have proved that the 
coefficients (a) satisfy the following recurrence relation [21 Lemma 
3.6]: 



(n,m) , 



2A; + 



a 



2n-l 



(n— l,m) 



a + 



f^/3f-~^V), (5) 



2m 



for k = 0, 1, m+n—1 and /3Q"''™''(a) = 0. From (jlj) and (jSD they derive 



(n,m) 



a) when < a < 1 and also that 



(n,m) I 



the positivity of 

for k < min(m,n). However, an explicit single sum formula for 



(n,m) 



a) is missing in their paper. 



Our main result is an explicit single sum formula for /3j^"'"*^(a), which 
provides actually the unique solution of the recurrence system (|5|) with 
the boundary condition (jlj). 



Theorem 1. For k = 0,1, ...,n + m, we have 



{n,m) I 



a) 



2{n+m—k) 

i=o 



-n+fc (l/2)n+m-fc(l/2)fc 
(l/2)„(l/2)„ 

n + m + 1 
2n + 2m-2k- j 



(6) 



— m + /c + j 
J 
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Equivalently, the coefficients can be written in terms of hy- 

pergeometric function: 

(i) ifk> \{n + m-l) /2] , then 

l3^-'-\a) = a'"^--'il - ^)-n+fe (V2Wfe(l/2). / n + m + 1 
> > (l/2)„(l/2)^ \2n + 2m-2k 

f -2m-2n + 2k, -m + k + 1 1\ . . 

""^^^ 2k-m~n + 2 'aj' 

(ii) ifk< [{n + m- 1)/2J , then 

x(-ir"-'.F,(-"-;-,;,-4-*;i). (8) 

We first derive some immediate consequences from tlie above tlieo- 
rem. Assume tliat n <m and < a < 1. 

(1) If < n, then tlie binomial coefficient { ^^^^} vanislies 
because n + m — 2A; — 1 > n — k — 1. Thus formula ([H]) implies 

that /3i"'"'V) = for A; < m A n. 

(2) If /c > n, then each term in ([B]) is a polynomial in a. Indeed, it 
is clear that (1 — a)"""'"'^ G Z[a], and 

• if > m, then 2n + m — k—j > 2n + m — k — 2{n + m — k) = 
k — m > for < j < 2{n + m — k), 

• if < m, then the binomial coefficient (^~'"+^+-'^ does not 
vanish only if —m + k+j < 0, that is, j < m — k, therefore, 
2n + m — k — j > 2n + m — k — {m — i) = 2n > 0. 

(3) From definition (|3]) we derive immediately the symmetry prop- 
erty: 

(3t"'\a) = pt''\^-^)- (9) 
This identity follows also from ([7]) (resp. ([8])) and Pfaff's trans- 
formation (see [U p. 94]): 

A B \ ,^ ,_A ^ f A, C - B 



with A = -2n-2m + 2k, B = -m + k + 1, C = -n-m + 2k + 2 
and X = 1/a if k > \{n + m — l)/2] (resp. A = —n — m — 1, 
B = n — k, C = n+m—2k and x = l/aiik < [(n + m — l)/2j). 

We shall prove Theorem 1 in Section 2. As applications of our Theo- 
rem 1, we shall derive a formula of Berg and Vignat [2] when m = n, a 
positivity result as well as a formula generalizing an integral evaluation 
of Boros and Moll [4] in Section 3. 

2. Proof of Theorem 1 
We shall first derive a double sum formula for the coefficients p^""'"^^ (a). 
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2.1. A double sum. Assume that m > n. Let 

-2n)kk\ 



(n) (-n)fc2^ 



By definition ([T]) we have 

m+n A; 

g„(an)g™((l-a)M) = n^' ^al'!^ (1 - a)'=--'\ (12) 

k=0 j=0 

Recall the connection formula due to Carlitz [S]: 

k 

= Y.5f\,{u), fc = 0,l,..., (13) 

where 

t^. { (fc+l)! (-l)>^-\2i)\ r ^- ^L. 

[ for < i < 

Plugging (IT^ into (IT^ and extracting the coefficient of qi{u) we obtain 

m+n k 

/5f ''"V) = E ^J'^ E «f - «)''^- (14) 

k=i j=0 

Replacing k hj k + i and j by n — j we obtain 

k=0 j=0 

(-f^ '^J)' (15) 

fc=o i=o ^ ^ 
^ (-n)n_j(-m)fc+i_„+ja"~-^'(l - a)^+^' 

{-2n)n-j{n - j)\{-2m)k+i-n+j{k + i - n + 

As fji'^'™'^ (a) = if z < n, replacing ihy n + m — i and assuming i < m 
we obtain, after some simplification, 

/'„N _ A-i (l/2)n.+m-t n+i/i _ \m-i_(^^^^_2tilL 

/^n+^-.(a)-4 (i/2)„(i/2)„« (1 «^ 



X 



^ (m + n - z + 2)A,.(-n - m + z - l)fc 

A;>0 



E(-i)* 



k\ 



y (-n),(-z)fe+,(n+ l),a--+^(l - l/a)^+^- 
By binomial formula, we have 



a-^+'^il - l/a)^-+^' = E ^ ^ , ^^* a--+^- 
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Let I = k + j and r = i + t — kwe can write 

(l/2)„(l/2),„ " 11 «' 1^ 2i 

xE°"'E , (IT) 

^ ^ (-m - «)z(l + m - ^ ^ 

r=0 />0 ^ ' 

^ {-l)k{-l)r-i+k{m + n- i + 2)k{-n -m + i-l)k 
k\{r-i + k)\{n-l + l)k{-n-l)k 

For 2 = 0, 1 and 2, the above formula yields immediately the following 
explicit expressions: 



(l/2)n(l/2), 



o(n,m) / N _ n+Ui _ „\m-l ( V2) (l/2)n+m-l 



n + m + l\ /n + m+l\/n\ fn + 1 

^ ' ' 2 j " V 1 / Vl/'' ^ V 2 



and 



= - a)--2(l/2)2(l/2)„+™_2 



(l/2)n(l/2) 



m 



^ + m + 1\ /n\ _ /n + m + 1\ fn + 1\ ^.3 



2 

n + 2' 



Further computation led us to conjecture the formula (jH]), which is 
equivalent to ([7]) and ([8]). 

Our proof consists of verifying that the conjectured formula ([6]) does 
satisfy the boundary condition (jl]) and the recurrence relation (Q. 



2.2. The boundary condition. We first show that the formula 



reduces to (jl]) when m = 0. Let cl"''(a) = Pl^'^\a) for < A; < n. 



Clearly we have ch\a) = a" from ([6]). Assume that < A; < n — 1, 
using the symmetry ch\a) = — a), we derive from ([6]) that 
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As (^^"^•') = if j > n — k, replacing j by n — A; — 1 — j we obtai 



am 



cr(a) = (l-a)a'^(l/2)„_. 



(l/2)„, 

(n-fc-l)Afc , I 1 \ / 1 ■ 

y (-ir^-^-Y V )(i-a) 



X 

Since 

(1/2), _ (^) / n + 1 \ _ fn + 1 



(l/2)„- 



^(i/2)„ g)' V^-A; + i+j7 V^-j 

and 



,n-fc-i-,^ -1-J A (j + l)(j + 2)...(n-fc-l) 



n — k — 1— jj [n — k — 1 — j] 

n — k — V 



-1 



J 

we derive the desired formula dl]). □ 

2.3. The recurrence relation. It remains to verify for the con- 
jectured formula ([6]). One of Gauss' relation for contiguous hypergeo- 
metric functions [3 p. 71] reads: 

C2Fi(^'^'^^ -z^ (18) 

This can also be verified straightforwardly by comparing the coefficients 
of z'' for each A; G N. Similarly we have the following contiguous 
relation: 

B{l-z),F,( "^^^^^^ ■,z) (19) 



{B - C)Z2F^ ( C '+ f ' ^) + ^2^1 ( 



A-1, B 
C 



Now, the recurrence relation ^ is satisfied by the conjectured formula 
([6]) since 

• for > \{n + m — l)/2] , this is equivalent to (ITSj) with z = 1/a 
and A = -2m - 2n + 2fc + 2, B = -m + k + 1, C = -n - 
m + 3 + 2k. 

• for k < \_{n + m — l)/2j , this is equivalent to (|T9l) with 2; = 1/a 
and A = —m — n, B = —n — k — 1, C = n — m — 2k. 
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As a byproduct, comparing and (fT7|l we derive immediately the 
following double sum identity. 

Corollary 2. For all positive integers i, m, n and r = 0, 1, . . . , 2i we 

have 

, min{i,n} mm{l,l+i—r} / ■\ / \ / -. \ 

m + i\ ^ ^ [-i)i[-n)i[n + l)i 



2i I ^ ^ (—m — i)i(l + m — i)il\ 

/=max{0,i— r } A;=max{0,i— r} 

{-l) k{-l)r-i+k{m + n-i + 2)k{-n -m + i -l)k 

X 



k\{r - i + k)\{n - I + l)fc(-n - l)k 
m + n+l\fn — i + r^ 
2i — r / \ r 



(20) 



It would be interesting to find a direct standard proof of fl2U]) using 
the hypergeometric function theory. 



3. Applications 

3.1. A positivity result. Theorem 1 implies immediately the follow- 
ing positivity result. 

Corollary 3. Suppose that n < m. Then the coefficient /3^"'™^(a) is 
positive for n < k < m. 

Proof. For (ra + m — l)/2 < < m, we have —m + + 1 < 0, —2m — 
2n + 2/c < and 2/c — m — n + 2>0. Hence each summand in ([7]) is 
positive and so is /3^"'™''(a). For n < k < (n + m — l)/2, we derive from 
(El) and ([TOD that 



= a'-{l - «)"^-''^'^^f7iyf^^ (21) 
n — k — 1 \ f —m — n — 1, m — k 1 

2^1 



n + m-2k-lj \ n + m -2k '1-a 

In this case, we have —m — n — 1<0, m — k>0 and m + n — 
2k > 0, therefore each summand in ( 1211) is positive and so is /3^'"^\a). 
Summarizing the above two cases, we obtain the desired positivity. □ 

3.2. A formula of Berg and Vignat. Let /3f"^(a) be the coefficients 
defined by 

n 

g„(au)g„((l - a)u) = ^ /3t\a)qn+i{u). (22) 

i=0 
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Corollary 4. For < i < n we have 

^ = (4«(l-«))-(-n).(n + p. 



Proof. By definition we have (i'f'\a) = /3^'^"''(a). It follows from 
that 



X 2i"l 



(l/2)„ V2^-2z 
-2n + 2ii, i + 1 1 
2z + 2 'a 

Applying the quadratic transformation formula [1, p. 127] : 



2 



2' 



we obtain 

A^")(a) = a^"-^(l-ar(l/2) 



W^.^_.2„-^/, /ox in + l/2)if2n+l 



(l/2)„ V2^-2z 



Reversing the order of summation in the last 2-^1 and using {x)n 

(-x-n+l)fc 



-1)^7^#TTT- we recover ([23]). □ 



3.3. An integral evaluation. Recall [2] the Student t-distribution 
with parameter u > 0: 



For z/ = n + i with n G N we say that the distribution has 2n + 1 
degree of freedom. The Fourier transform of the density fy is an even 
function, which for ^ > 0, can be expressed as 

/oo nl — v 

where K^, is the modified Bessel function of the third kind, also called 
the Macdonald function. 

For 1/ = + I it is known that k^{u) = e~"g„(u) for m > 0. Thus 
equation is equivalent to 

k=nAm 
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where < a < 1 and * is the ordinary convolution of densities. 

For a = replacing x by | and multiplying by ^ on both sides, we 
obtain 

n+m ^ ^ 



2 " ""^5^2 

k=nAm 



I 5 



which is equivalent to the following integral evaluation. 

, ^ m+n A n{n,m)n- 

dy - Yl 



(1 + y^r+'il + ix- 2 ^^^^ (1 + 



By Theorem 1 with a = 1/2, the above formula can be simplified as 
follows. 

Corollary 5. For positive integers m and n we have 

n!m!2"+" dy _ n 7^"'"^^ 



where 



2{n+m—i) 1 \ / • ■ 

Tk = 2^ (-2)' ' ' ' 

j=0 



2n + 2m — 2i — j J \ j 



When m = n, as shown by Berg and Vignat [3], the above formula 
yields a new proof of an integral formula of Boros and Moll |1] . 
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